Dynamical Behavior of Spins in the Rare-Earth Kagome Pr3Ga5SiOi4 
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We report on the use of ^^'^'^Ga nuclear magnetic resonance to probe spin dynamics in the rare- 
earth kagome system Pr3Ga5SiOi4. We find that the spin-lattice relaxation rate ^^l/Ti exhibits a 
maximum around 30 K, below which the Pr'^'^ spin correlation time r shows novel field-dependent 
behavior consistent with a field-dependent gap in the excitation spectrum. The spin-spin relaxation 
rate ^^l/Tg exhibits a peak at a lower temperature (10 K) below which field-dependent power- 
law behavior close to is observed. These results point to field-induced formation of nanoscale 
magnetic clusters consistent with recent neutron scattering measurements. 

PACS numbers: 75.10.Hk, 75.30.Gw, 76.60.-k 
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Various two-dimensional (2D) triangular lattices, 
including kagome systems, exhibit intriguing low- 
temperature spin dynamics induced by a geometrical 
frustra,tiOT [Jf^i, i, i, i, i, 0, i, 1, M 
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The prototypical kagome examples 
are ZnCu3(OH)6Cli2 [Cu^+ (S=l/2)], SrCr8_,Ga4+.Oi9 
[Cr3+ (S=3/2)], and the jarosite (D30)Fe3 (804)2(00)6 
[Fe3+ (S=5/2)] i, i, i]. These are all 3d transition 
metal based compounds in which the exchange interac- 
tions and geometrical frustration are of dominant impor- 
tance while Dzyaloshinsky-Moriya interactions, single- 
ion anisotropics, and ofF-stoichiometry serve as small per- 
turbations. Nonetheless, investigation of the predicted 
spin-liquid state is complicated by the presence of these 
perturbations. 

For rare-earth (RE) compounds, the situation might be 
less complicated: the single-ion anisotropics govern the 
overall magnetic behavior, while frustration-induced spin 
dynamics is expected to emerge at energies well below the 
crystal field splitting. The recently discovered RE based 
kagome compounds R3Ga5SiOi4 (R=Nd or Pr) do show a 
disordered state at low temperatures [l7 , H, 1^, 2^ ■ The 
cooperative magnetic correlations depend on the RE ion 
involved and the crystal field environment. However, this 
issue has not yet been fully addressed. It is therefore de- 
sirable to study the spin dynamics of Pr3Ga,5SiOi4 and to 
compare the behavior with its isostructural counterpart 
Nd3Ga5SiOi4 which has different magneto-crystalline 
anisotropy poj . 

Pr3Ga5SiOi4, which belongs to the langasite family, 
has a trigonal crystal structure (P321) with lattice pa- 
rameters a = 8.0661(2)1 and c = 5.0620(2)1. The RE 
Pr'^+ ions (5/^ : J = 4) are networked by corner sharing 
triangles in well-separated planes to form a distorted lat- 
tice (see the inset of Fig. [1]) topologically equivalent to the 
kagome lattice. Analysis of the specific heat suggests that 
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FIG. 1: (color online). NMR field-scan spectra of 
Pr3Ga5SiOi4 at 92 MHz showing quadrupolar split ®^Ga 
(J = 3/2) and ^^Ga (7 = 3/2) components for the three non- 
equivalent Ga sites. Resolved peaks are labelled central (C) 
or satellite (S) components. The left inset shows the enlarged 
version of the portion of the graph in box. The right inset de- 
picts the crystal structure showing three kagome planes. Two 
Ga sites lie between these planes and the third site is in plane 
(see Ref. 



the crystal field energy level structure of Pr'^^ consists of 
three singlets with gaps Ai = 25 K and A2 = 68 K, and 
A3 ~ 780 K At high temperatures the magnetic 

susceptibility x obeys the Curie- Weiss law with a Curie- 
Weiss temperature of ^cw = —2.3 K that points to weak 
antiferromagnetic (AF) correlations between the Pr'^+ 
spins [ll]. At zero external field Pr3Ga5SiOi4 shows (i) 
no long-range magnetic ordering at temperatures down 
to 35 mK, (ii) a dependence of the specific heat at low 
T, and (iii) spin excitations consistent with a highly de- 
generate state [13. These suggest that Pr3Ga5SiOi4 has 
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a spin-liquid-like ground state corresponding to a high 
frustration index f = ^cw/Tc > 66. Application of an 
external field parallel to the c-axis leads to the forma- 
tion of nanoscalc magnetic clusters, whose size increases 
with H. This behavior is accompanied by reduction of the 
component of the specific heat and the opening of a 
spin gap in the excitation spectrum [Tst - 

^SR and NMR measurements on the sister RE kagome 
compound Nd3Ga5SiOi4 reveal that a disordered state 
persists down to low T in zero field [li,[2^. Spin fluc- 
tuations are suppressed by an applied field and for H = 
0.5 T a field-induced transition is found at 60 mK. Al- 
though the two isostructural compounds R3Ga5SiOi4 (R 
= Pr'^+ or Nd'^"'") share, to some extent, common physics, 
detailed cooperative spin dynamics varies with the RE 
ions. This is because Pr'^+ and Nd'^"'" have different crys- 
tal field splittings, single-ion anisotropics, and exchange 
interactions leading to a substantial difference in low- 
energy spin dynamics both at zero field and in external 
field. 

In this Letter we present ^^'^^Ga NMR lincshapc and 
relaxation rate measurements made on Pr3Ga5SiOi4 as 
a function of temperature in applied fields in the range 
2 — 17 T. The objective is to follow the spin dynamical be- 
havior with temperature in the spin liquid regime and to 
determine how the dynamical behavior depends on large 
applied magnetic fields. We find that spin-spin relaxation 
rate is sensitive to field-induced short-range magnetic or- 
dering. This enables us to separate frustration-induced 
cooperative phenomena from single-ion physics. 

Single crystal samples of Pr3Ga5SiOi4 were grown us- 
ing the traveling floating-zone technique and character- 
ized by X-ray diffraction as described previously 
The NMR spectra and relaxation rate measurements 
were obtained using a pulsed spectrometer with quadra- 
ture detection. 

Figure [1] shows the NMR field sweep spectrum ob- 
tained by integrating spin-echo signals. Multiple peaks 
are found corresponding to three non- equivalent sites 
for the / = 3/2 ^^Ga and ''^Ga isotopes each hav- 
ing quadrupolar splittings which give rise to a central 
line and two satellites due to non-cubic site symmetry. 
The resulting spectrum consisting of eighteen overlap- 
ping hues is similar to that of Nd3Ga5SiOi4 Two 
Ga sites lie between the kagome planes and site 3 is ran- 
domly occupied by Ga'^"'" and Si^"*" ions [23|. The princi- 
pal components of the spectrum arc denoted as C(l) and 
C(2) for each isotope. 

Figure[2ta) plots the measured spectral NMR shift ^'^K 
for _ff = 9 T along c versus x\\ with T as the implicit pa- 
rameter while Fig. [Ifb) shows the spectra in a stacked 
plot. The linewidth increases significantly as T is low- 
ered below 100 K reaching a plateau value of 0.22 T below 
10 K. For T > 30 K a linear relationship between ^^/i 
and x\\ is found. However, departures from this relation- 
ship become apparent below 30 K where spin correlation 
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FIG. 2: (color online), (a) ®^Ga NMR shifts measured m an 
applied field of 9 T along the crystal c-axis plotted versus 
the magnetic susceptibility with temperature as the implicit 
parameter. For T < 30 K the plot shows departures from a 
simple linear relationship, (b) "^Ga field-swept NMR spectra 
at 95 MHz as a function of T. 



effects become increasingly important. 

In order to study the Ga'^+ ion spin dynamics, spin- 
lattice (l/Ti) and spin-spin [I/T2) relaxation rates were 
measured as functions of temperature for the central 
^^Ga spectral component corresponding to ^^C(l) (see 
Fig. El^b)). Relaxation rate measurements were made 
in several different applied magnetic fields directed par- 
allel to the c-axis and the results are given in Figs. [3] 
and m Similar results were found for the other spec- 
tral components. At low T the magnetization recovery 
curve showed stretched exponential form M — Mq[\ — 
exp(— t/Ti)''] where (3 w 0.7 indicating a distribution 
of relaxation rates. In contrast to the sister compound 
Nd3Ga5SiOi4 [i^l and other transition metal based frus- 
trated 2D antiferromagnets such as NiGa2S4 
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which wipe-out of the NMR signal occurs at low temper- 
atures due to low frequency magnetic fluctuations, the 
NMR signal in Pr3Ga5SiOi4 could be observed in all ap- 
plied fields over the temperature range 300 mK to 290 K. 
As a result, the relaxation rates could be measured over 
the entire temperature range showing well defined max- 
ima below 30 K, as seen in Figs. [3] and |4l permitting 
detailed analysis of the data. 

As noted above, magnetic contributions to the specific 
heat are accounted for by assuming low-lying crystal field 
split states for the Pr'^^ ions Following this model, 
we may expect that, for T < 30 K, the temperature de- 
pendence of the correlation time due to transitions be- 
tween the ground state and the first excited state should 
be given by r = tqc^^^'^ with Ai the lowest energy gap 
and the pre-exponential factor tq ^ 10^^^ s. This pre- 
dicts that for T < Ai the correlation time for spin fluctu- 
ations will increase rapidly with decreasing T. Behavior 
of this kind has been found in /iSR and Ga NQR relax- 
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FIG. 3: (color online). '^^Ga 1/Ti and I/T2 for PraGasSiOw 
as a function of T at 9 T, together with magnetic specific 
heat at 9 T whose peak is close to the maximum of ^^1/T2. 
The similarity in behavior of ®^ 1 /T2 and Cmag is striking and 
points to a common underlying mechanism. Inset: log-log 
plot of ^^1/T2 at different fields. Notice that the broad max- 
imum sharpens as the field is increased and below the peak 
the behavior is close to shown as the dotted line. 



ation rate measurements in the quasi-2D AF NiGa2S4 
14| . However, the present relaxation rate behavior can- 
not be accounted for using the Arrhenius expression with 
a field-independent energy gap. The slope of the low T 
region of the log- log plot of 1/Ti versus T in Fig. [3] (in- 
set) gives a field-dependent gap AyvA/i? which is discussed 
below. 

Spin-lattice relaxation is attributed to fluctuating hy- 
perfine fields, produced by the electron moments on 
nearby Pr'^+ ions, which induce nuclear spin state transi- 
tions. The fluctuating local field H]^ at a Ga nuclear site 
is mainly due to dipolar interactions with neighboring 
electron spins and any transferred hyperfine interaction 
plays a subsidiary role. The dipolar interaction induces 
nuclear transitions at the Larmor frequency loj while 
an isotropic hyperfine interaction will induce mutual 
electron-nucleus transitions at the much higher frequency 
|a;/±ws| where tus is the electron frequency. For long cor- 
relation times the dipolar process is dominant and we can 
neglect any contributions due to the hyperfine coupling. 
Assuming that the spin correlation function decays ex- 
ponentially, the spin-lattice relaxation rate may be writ- 
ten as l/Ti = i+Ii|t^ ) where T2 is the transverse 
correlation time for electron spins and C± = lliH"]) 
is the transverse component of the local field [2l|, l22j . 
Similarly, we have \/T2 = \C! i-ij^^jr;^) + C'||Ti with ri 
being the longitudinal correlation time, C\\ = jj{H'u), 



and the z-component of the local field. The expres- 
sion for 1/T2 holds provided ti < 1/Alo where Auj is 
the NMR linewidth of the selected spectral component. 
For Ti < l/Aw, I/T2 is expected to saturate. The in- 



troduction of both transverse and longitudinal electron 
spin correlation times allows for the possibility of differ- 
ent mechanisms being of dominant importance for ti and 
T2, respectively. A maximum in l/Ti occurs for W/T2 = 1 
while for higher T, in the short correlation time case 
(cjir2 < 1) the dipolar mechanism leads to 1/Ti ~ C±T2 
and I/T2 ~ ^C±T2 + C||Ti. Inspection of Fig. [3] shows 
that for T >" 30 K, the 1/Ti and I/T2 curves lie close 
together suggesting that the transverse fluctuations are 
of dominant importance, corresponding to C± > Cy , for 
both relaxation rates in this interval. It is likely that 
Ti — T2 at high T. Use of the condition for the maximum 
in 1/Ti permits T2 values to be obtained as a function 
of T as shown in Fig. 2) Assuming the Arrhenius rela- 
tion holds in the high-T (T > 30 K) region the fitted 
curve gives tq ~ 10^^^ s and the energy gap for spin ex- 
citations A = 98 K. This is consistent with the crystal 
field splittings (Ai + A2) given by the specific heat re- 
sults [18[. The energy gap obtained from the slopes of 
the curves in the low T region of Fig.[51 denoted A^mr, 
is clearly field-dependent with values plotted versus H in 
the upper inset. The fitted curve in this plot has the form 
An MR = Ao + aH where the slope a ^ g^B with /is the 
Bohr magneton and g = 3.32 close to the g value for the 
Pr'^+ ion. The magnitude of the zero-field gap, Aq = 3.5 
K (see Fig. |4] caption) , obtained from the low tempera- 
ture (T < 10 K) NMR relaxation data is much smaller 
than the high-T NMR value or that from the specific 
heat results. We note that integrated inelastic neutron 
scattering data at 35 mK [3l give a field - dependent 
spin gap in the excitation spectrum strikingly similar to 
the gap derived from the NMR correlation time behav- 
ior as shown in the upper inset in Fig. 2) It is likely that 
field-suppressed magnetic fluctuations are responsible for 
the observed fleld-dependence. 

While 1/Ti decreases at temperatures below the max- 
imum shown in Fig. [31 1 /T2 continues to increase as T is 
lowered before passing through a maximum and then de- 
creasing dramatically below 7 K. The behavior is strongly 
field-dependent and this again points to field suppression 
of magnetic fluctuations. Fig. [3] compares the behav- 
ior of 1 /T2 at 9 T with that of the specific heat at 9 T 
with the scales adjusted to allow comparison. The mag- 
netic speciflc heat Gmag{T), which is obtained from the 
measured Cp{T) by subtracting the lattice contribution 
using La3Ga5SiOi4 as a reference, shows behavior for 
H = 0.15 T. The entropy saturates at i?ln3 as expected 
for a system with three low- lying crystal fleld states. The 
temperature dependence of l/r2 is strikingly similar to 
that of Cmag (T) as seen in Fig. [3] (for H = 9 T) and this 
similarity in form points to a common underlying process 
involved in both observed behaviors linked to the spins 
entering the singlet ground state at low T. Similar be- 
havior is observed in NiGa2S4 where the estimated spin 
"freezing" temperature is coincident with the peak in 
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but NMR wipeout prevents detailed comparison 
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FIG. 4: (color online). Temperature dependence of the trans- 
verse spin correlation time T2 at different fields extracted 
from vs T plot (see lower inset). At high tempera- 

tures, the gap A ~ 98 K obtained from the slope is field- 
independent. Below 30 K, the gap Anmr has field depen- 
dence (see upper inset) where the dashed line corresponds 
to the fit Anmr = Aq -)- aH with a = g^B = 3.32/^_b and 
Aq — 3.5 K. The field-dependence of the spin gap in the ex- 
citation spectrum as derived from 35 mK inelastic neutron 
scattering results (Ref. [T^ ) is shown for comparison. 



('--^ 6 to 7 in-plane lattice spacings) for = 9 T [18j . 
The present relaxation rate results suggest that at low 
temperatures the spin correlation time r becomes long 
as shown in Fig. 3] and a possible reduction in the c-axis 
component of the dipolar field at Ga sites results from 
in-plane short range ordering leading to the anomalous 
behavior of 1/T2. 

In conclusion, we have investigated the spin dynam- 
ics in the rare-earth kagome system Pr3Ga5SiOi4 via 
^^Ga NMR measurements over a complete temperature 
range without the NMR wipcout effect seen in the sis- 
ter compound Nd3Ga5SiOi4. The correlation time r for 
spin fluctuations extracted from the spin-lattice relax- 
ation rate values exhibits novel features; r increases with 
decreasing T and below 30 K the results arc consistent 
with a field-dependent energy gap in the excitation spec- 
trum. The spin-spin relaxation rate shows a maximum 
close to that in the specific heat and the form of the 
temperature dependence of these two quantities is very 
similar below 10 K. This result warrants theoretical at- 
tention. The drop in 1/T2 below 7 K is attributed to 
a decrease in the local field at Ga sites linked to field- 
induced nanoscale ordering of the electron spins in the 
a-b plane. 
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tive Agreement No. DMR-0084173, EIEG grant, by the 
State of Florida, and by the DOE. 



with I/T2 [14|. Pr3Ga5SiOi4, however, allows this com- 
parison to be made over a complete temperature range. 
Below 7 K occupation of the singlet ground state by the 
electron spins rapidly increases with decreasing T and 
both I/T2 and C'mag drop to low values. The behavior of 
I/T2 below 7 K may be accounted for by decreases in ti 
and/or While it is unlikely that ri will shorten at low 
T, local ordering of spins could result in a reduction in 
leading to a change in the "static" field contribution 
to I/T2. This effect will be particularly marked if spin 
ordering occurs in the a-b plane. The NMR spectra for 
T < 10 K show little change in linewidth and the Knight 
shift is roughly constant (Fig. [5]). 

The quadratic T-dependence of Cmag for = T is 
interpreted as evidence for gapless Goldstone modes. For 
H > Q T, Cmag{T) has a minimum at 80 mK with the 
low- T upturn ascribed to nuclear contributions to the 
specific heat [l^. Elastic neutron scattering measure- 
ments for T < 1 K show that nanoscale ordering occurs 
in the presence of an applied field consistent with 2D 
short-range order and give a correlation length of 29 A 
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